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ABSTRACT 



Context. Diffuse envelopes around Mira variables are among the most important sources influencing the chemical evolution of galax- 
ies. However they represent an observational challenge because of their complex spectral features and their rapid temporal variability. 
Aims. We aimed at constraining the exact brightness distribution of the Mira star T Lep with a model-independent analysis. 
Methods. We obtained single-epoch interferometric observations with a dataset continuous in the spectral domain (A = 1.5 - 2.4yum) 
and in the spatial domain (interferometric baselines ranging from 1 1 to 96 m). We performed a model independent image reconstruc- 
tion for each spectral bin using the MIRA software. We completed the analysis by modeling the data with a simple star+layer model 
inspired from the images. 

Results. Reconstructed images confirm the general picture of a central star partially obscured by the surrounding molecular shell of 
changing opacity. At 1.7 yum, the shell becomes optically thin, with corresponding emission appearing as a ring circling the star. This 
is the first direct evidence of the spherical morphology of the molecular shell. Model fitting confirmed a spherical layer of constant 
size and changing opacity over the wavelengths. Rough modeling points to a continuum opacity within the shell, in addition to the 
CO and H2O features. Accordingly, it appeared impossible to model the data by a photosphere alone in any of the spectral bins. 

Key words, techniques: interferometric - stars: AGB and post-AGB - stars: atmospheres - stars: individual: T Lep - stars: mass-loss 



1. Introduction 

Among the Asymptotic Giant Branch members, Mira stars are 
low-mass (1M ), large-amplitude (AV ~ 9), long -period vari- 
ables (« lyr), evolving toward the planetary nebula and white 
dwarf phases. Their important mass-loss rate, of the order of 
10 _6 M o /year, significantly affects their evolution and is one 
of the most important sources for the chemical enrichment of 
the interstellar medium. A better understanding of these late 
stages of stellar evolution may help shed light on this impor- 
tant player in the chemical evolution of galaxies. Due to their 
large diameters and high luminosities, Mira variables are a fa- 
vorite target for observations at high angular resolution in the 
optical and near-infrared (NIR) wavelengths. The wealth of 
aperture-masking and long baseline interferometric information 
has strongly advanced the study of molecule and dust formation 
dPerrin etalJll999t [T hompso n et al.ll2002t IWoodruff et al.l[2004t 
IOhnakal2004l:IIreland & Scholzll2006l:lRagland et al.ll20061) . The 
global picture is an onion-like structure with molecular layers 
and dust shells surrounding the photosphere, the latter animated 
by various activities (pulsations, convection, wave-shocks) cru- 
cial to explain the mass-loss process and the structural evolution. 

To fully characterize the structure of a Mira star's atmo- 
sphere, one would need the complete intensity map at all wave- 
lengths and all pulsation phases, recorded over a large number 
of pulsation cycles. The majority of recent studies concentrated 
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on the spectral and/or the phase dependence. They overcame the 
issue of the exact brightness distribut ion by using either inte- 
grated quantities (su ch as diameter, Mill an-Gabet e t al. 2005; 
Woodruff et al. 2008), or the direct com parison of a few interfer- 
ometric measurem ents with geometric dMennesson et alj|2002l : 
Perri n et alj|2004l) or dynamic models of pulsatin g atmospheres 
dOhnaka et al.ll2006t IWittkowski et alj|2007l 120081) . To study the 
brightness distribution, the difficulty lies in collecting a sufficient 
amount of spatial information in a period of time shorter than the 
typical life-time of the expected structures, and within spectral 
bins small enough to spectrally resolve the molecular bands (for 
which completely different intensity maps are expected). 

In this work, we present a multi-wavelength, spatially re- 
solved observation of the single star T Lep. This Mira variable 
has a spectral type of M6e-M9e and V -band magnitude vary- 
ing from 7.3 to 14.3 dSamus et al.l 2004). For the ephemeris, we 
use a Modified Julian Day of last maximum brightness 7q = 
54446 days and a pulsation period of 380 days (based on the 
last cycles observed by AAVSO), instead of the 368 days from 
Whit elock et al.l d2000l). The mass- l oss ha s been estimated to be 
7.3 x 10~ 7 M o /year by lLoup et all d!993l) . The most prominent 
spectral features expected in the NIR comes from water vapor 
and also CO. The shapes and widths of these molecular bands 
depend on the stellar phase. 

We present the dataset in Sec. |2l with a special emphasis 
on the observations and data reduction. Image reconstruction is 
presented and discussed in Sec. [3] We mod el the data in Sec. [4] 
with a simple star+layer model inspired bv lPerrin et alj d2004l) . 
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Fig. 1. Example of observation sequence, on the baseline G0-H0 
and for the spectral bin 1 .59 /urn. The interpolated TF is repre- 
sented. 

Table 1. Observation log. 



Date 


T Lep phase 


Baselines 


Spectral range 


2008-10-10 


0.80 


G2-G0-K0 


1.5087 -2.3191yum 


2008-11-04 


0.86 


D0-H0-G1 


1.4287 - 2.44 13 fim 


2008-11-06 


0.87 


A0-D0-H0 


1.4287 -2.423 8 yum 


2008-11-10 


0.88 


E0-G0-H0 


1.4287 -2.441 3 yum 



2. Observations and data reduction 

Data w ere collected at the Ver y Large Telescope Interferometer 
(VLTI. lHaguenauer et alj|2008[) with the spectrograph AMBER 
dPetrov et al.l T2007F covering simultaneously the J-, H- and K- 
bands with a spectral resolution of R ^35. Even though J-band 
fringes have been properly recorded in several observations, we 
decided to discard them from the analysis of this paper since the 
data quality is significantly worse than for longer wavelengths. 
The majority of the o bservations have been obtained by the 
FINITO fringe-tracker dLe Bouquin et a l. 2008). When the con- 
ditions were unstable, we used FINITO in group-tracking mode 
(instead of phase-tracking). This mode reduces the sensitivity to 
unstable atmospheres at the cost of a clearly reduced instrumen- 
tal contrast. We used 4 configurations of 3 Auxiliary Telescopes 
(ATs) each. Data for the A0-D0-H0, D0-H0-G1 and E0-G0-H0 
configurations were obtained within a few days of each other 
while data with G2-G0-K0, extracted from the archive, was ob- 
tained about 1 month before as shown in Table [1] 

Concerning the wavelength tables, we performed a rough 
absolute calibration by cross-correlating our observed spectra 
with an atmospheric model. We found a systematic offset of 
0.08 + 0.02 ^m with respect to the initial AMBER table. 

Observations of calibrator stars with smaller angular diame- 
ters were interleaved with observations of T Lep. We were able 
to reach an average observation frequency of about 25min per 
calibrated point. This appeared to be crucial to properly sam- 
ple the instrumental and atmospheric transfer function (TF). 
This also allowed us to manually discard the obvious out- 
lier measurements without degrading the spatial coverage of 
our dataset too much. Additionally, we discarded all observa- 
tions where the FINITO locking ratio was lower than 60%. 
Raw visibility and closure phase values were computed using 
the latest public version of the am dlib package (version 2.2, 
iTatulli & A MBER consortium! 120071) and the yorick interface 
provided by the Jean-Marie Mariotti Center. We used our own 
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Fig. 2. Projected baselines (UV-plane) in meters of the complete 
AMBER dataset on T Lep. 



software to calibrate the instrumental and atmospheric TF (see 
Fig. Q] for an example). We estimated the uncertainty on the TF 
with the dispersion between the consecutive points obtained for 
each observation of a calibration star. This dispersion generally 
dominated the other sources of error (mainly the diameter uncer- 
tainty of calibrators) and was propagated to the uncertainties on 
the calibrated visibilities and phases obtained for T Lep. Fig. [2] 
shows the final UV-plane coverage of all observations that suc- 
cessfully passed all steps of the data reduction and calibration 
quality control. The east-west direction is favored because of the 
geometry of the telescope triplets used. Visibility curves are dis- 
played in Fig. [3] 

Looking at the interferometric data, we found that the gen- 
eral picture of T Lep significantly departs from a simple disk or 
Gaussian. It is composed of at least two different features: one 
has a characteristic size 5 - 7mas and the other 11 - 16mas. 
The exact size of these components and the flux ratio depend 
on the wavelength. This is qualitatively in agreement with a 
complex stellar atmosphere surrounded by molecular shells with 
wavelength-dependent opacity. Additionally, the non-zero clo- 
sure phase (departure from 0-180 deg up to 90 deg in the H-band 
and 30 deg in the K-band) indicated the presence of asymme- 
tries either in the shells or in the photosphere itself. Interestingly, 
we were not able to reproduce correctly the interferometric data 
with simple geometrical models such as Gaussian + Gaussian, 
or limb-darkened disk (LDD) + Gaussian, even when disregard- 
ing the closure phases This was our main motivation to attempt 
to reconstruct an image. 



3. Image reconstruction 

Our dataset is perfectly suited to spectro-spatial inversion, where 
the spectral and the spatial dimensions are treated together as- 
suming a certain number of degrees of freedom in each space. 
Unfortunately the necessary algorithms are not yet available, 
and we had to solve the problem independently for each spec- 
tral bin. We performed the image reconstru ction with the MIRA 
software (Thiebaut 20081; ICotton et alj|2008l) . In spite of the ap- 
parent completeness of the dataset, we expected that it would 
not be straightforward to produce images: (i) the geometry of 
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Fig. 3. Visibility profiles V(Bi) from AMBER/VLTI, compared to the star+layer model of iPerrin et al ] (12004 . The horizontal axes 
are the spatial frequencies B^, marked in meter per micron, the vertical axes are the linear visibilities. 



the VLTI array makes the east-west direction privileged, (ii) the 
necessary use of closure phases makes the phase information 
three times less rich than that of the visibilities. Therefore, a 2- 
step image reconstruction strategy was specifically elaborated. 
The first step consisted of building a radially symmetrical image 
from the data. Only the baseline length was considered. The am- 
plitude and sign of these pseudo data (which are all real due to 
the assumed symmetry) were derived from the measured power 
spectrum and sign of the closure phases respectively. We then 
used MIRA with a strong smoothness regularization to recon- 
struct the first stage images. The second step uses these bright- 
ness distributions as a support for a quadrati c regularization of a 
2D brightness map (Le Besnerais et al. 2008). Each of the spec- 
tral channels was processed separately. Resulting final images 
are presented in Fig. |4] 

The recontructed images clearly highlight two components. 
The brighter central component corresponds to the photosphere, 
with a linear ra dius of 100 R e (assuming a parallax of 5.95 mas 
from Lee 2007). The second component is the molecular layer 
(molsphere), characterized by the high spectral dependency of 
its apparent morphology. Within the molecular absorption bands 
(1.5 - 1.6^m, 1.75 - 2.15/mi, 2.3 - 2.5 /-im), the shell is opti- 
cally thick and appears, projected toward the observer, as a cir- 
cular disk. Between these bands, absorption by molecules is not 
as strong, and the shell becomes optically thin. It appears as a 
projected ring encircling the photosphere (1.6 — 1.75 /im). The 
lack of data points in the second lobe at 2.2 /mi (see Fig. |3) pre- 
vented the algorithm from reconstructing the inner gap between 
the molecular shell and the photosphere. Instead, it resulted in a 
decreased contrast between the molsphere and the photosphere. 

Each spectral bin has been reduced and imaged indepen- 
dently. These images are in good agreement. It provides a strong 
indication that the MIRA algorithm is not inverting the statistical 
noise present in the data. Possible correlations may still originate 
from the calibration noise (atmospheric turbulence) and the use 
of common calibrator stars. However, these sources of errors are 
unlikely to give rise to such a coherent pattern as observed in 
Fig. 2] These images are a strong argument for a shell-like ge- 
ometry of the molecules within the atmosphere. 



4. Intensity profile modeling 

In addition to reconstructed images of the molecul ar shells, we 
used a simple radiative transfer model developed bv lPerrin et al.l 
(2004) to derive the physical parameters of the system. This 
model consists of a uniform disk star that emits radiation with 
a black body spectrum (O*, T+), surrounded by a spherical layer 
with no geometrical thickness. The layer is characterized by a 




X (um) 

Fig. 5. Size and optical depth of the layer surrounding the pho- 
tosphere, as extracted by a fit of the visibility curves (Sec. 0J. 
The dashed curve is a fit with H2O and CO opacities assuming 
Ti = 1800K (T. Verhoelst, private communication). A contin- 
uum opacity of t ~ 0.2 was added. 



size, a temperature and an optical depth (O/, Ti, tj). Such an ana- 
lytic model can be compared to our wealth of observations with- 
out requiring time-consuming computations. We have also pre- 
ferred to use an ad-hoc geometrical model rather than a self- 
consistent hydrodynamical model in order to easily give some 
geometrical flexibility to the fitting process. Other types of mod- 
els will be discussed in further publications. We allowed the 
layer size and optical depth to be wavelength dependent, but 
not its temperature nor the stellar parameters. We were able to 
reproduce the visibility curves with the following parameters 
= 3500 K, T, = 1800K, <D* = 5.8 mas and the chromatic 
quantities of Fig. [5] According to the good overlay of the vis- 
ibility curves in Fig. [3] all spectral features observed in Fig. 
are real. However, our fit shows important correlations between 
T+, T\ and the average level of t; Additional spectroscopic 
inputs may be needed to provide definite measurements of these 
quantities. 

The model is able to reproduce the wavelength dependency 
of the visibility break around B^ = 20. Our data strongly suggest 
that the star+layer geometry is a good proxy for the real inten- 
sity profile and strengthens the result we had obtained with our 
imaging technique. It also gives a strong argument suggesting 
that the opacity effect explains most of the apparent size changes, 
without requiring important changes in the layer size (which was 
wavelength-independent in previous studies). 
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Fig. 4. Reconstructed images of T Lep with the MIRA software for several AMBER spectral bins across the H and K bands. The 
interferometric beam size is displayed in the bottom-right part of each image. Spatial scale is in mas. The white circles in the first 
column represent the average radius for the molecular layer (<£/ ~ 15mas) and for the central star (O/ ~ 5.8mas), extracte d from th e 
modeling of Sec. [4] It corresponds to the respective diameters of 2.5 and 1 AU (assuming 5.95 + 0.70mas parallax from lLed l2007). 
The mean surface brightness ratio between the photosphere and the molecular environment is around 10%. 



Following computations of lVerhoelst et alj d2006l) . we mod- 
eled the spectral dependency of the optical depth as a sum of wa- 
ter vapor (mainly contributing in the range 1 .75 - 2.15 fan), and 
CO (below 1.75 /mi and above 2.3 /jm). However, we had to add 
a continuum opacity of constant value *0.2 to be able to fairly 
reproduce the observed curve. Accordingly, the optical depth of 
the layer is always greater than 0.4 over the complete spectral 
range : it is impossible to see the photosphere alone. The pres- 
ence of this continuum opacity, as well as the discrepencies in 
the fit, will require adapted modeling before being understood. 



The simple star+layer model was able to reproduce almost 
all the visibility points with V > 15% but clearly failed at higher 
spatial frequency. An important feature in the visibility curves is 
the increase of contrast with baseline observed around B,j = 45 
for the spectral bins around 1 .8 fim. It could be tempting to iden- 
tify it as a "visibility second lobe" following a "visibility null" 
(confirmed by the closure phases). This visibility lobe is appar- 
ently no longer observed around 1.9 /mi and above. The only 
way to continuously remove a visibility null is i) to fill it with 
unresolved flux or ii) to displace it to infinite spatial frequen- 
cies by smoothing the corresponding edge, and iii) to add im- 
portant asymmetries so that the visibility becomes a complex 
number around the null. The two first cases are not plausible in 
our case because the visibility around = 45 never reaches 
in the entire dataset. We conclude that the visibility profiles 
cannot be explained by an axisymmetric distribution continu- 
ous in wavelength. Undoubtedly, such an effect would have been 
missed with a less complete spectral and/or spatial coverage. 



5. Conclusions 

We believe this letter presents the first interferometric observa- 
tions of a Mira star with a dataset almost continuous in the spec- 
tral and in the spatial domain. Obtaining it in a short amount of 
time has been made possible by the use of the AMBER spectro- 
graph associated with the relocatable ATs of VLTI. 

Individual spectral channel image reconstruction was 
achieved using the MIRA software. It provides confirmation of 
the existence of a spherical molecular shell at a 1 AU of the 
photosphere. Between 1.5 — 2.4 //m we did not find any spectral 
region where the brightness profile of T Lep can be explained 
by a single Gaussian or a single LDD. Conversely, when using a 
simple radiative transfer model we never found an optical depth 
less than 0.4 for the layer surrounding the photosphere. Rough 
modeling points to a continuum emission within the molecular 
shell, in addition to the CO and H2O features. 
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